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ABSTRACT 


The origin and depletion of iron differ from all other abundant refractory elements that make up the 
composition of the interstellar dust. Iron is primarily synthesized in Type Ia supernovae (SNe Ia) and 
in core collapse supernovae (CCSN), and is present in the outflows from AGB stars. Only the latter 
two are observed to be sources of interstellar dust, since searches for dust in SN Ia have provided 
strong evidence for the absence of any significant mass of dust in their ejecta. Consequently, more 
than 65% of the iron is injected into the ISM in gaseous form. Yet, ultraviolet and X-ray observations 
along many lines of sight in the ISM show that iron is severely depleted in the gas phase compared 
to expected solar abundances. The missing iron, comprising about 90% of the total, is believed to 
be locked up in interstellar dust. This suggests that most of the missing iron must have precipitated 
from the ISM gas by cold accretion onto preexisting silicate, carbon, or composite grains. Iron is thus 
the only element that requires most of its growth to occur outside the traditional stellar condensation 
sources. This is a robust statement that does not depend on our evolving understanding of the dust 
destruction efficiency in the ISM. Reconciling the physical, optical, and chemical properties of such 
composite grains with their many observational manifestations is a major challenge for understanding 
the nature and origin of interstellar dust. 
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1. INTRODUCTION 


Interstellar dust condenses in quiescent stellar outflows 
during the asymptotic giant branch (AGB) phase of their 
evolution, or in the explosive ejecta of CCSNe. Its pres- 
ence in the ISM is manifested primarily by the extinc- 
tion of starlight, the diffuse infrared (IR) and millime- 
ter emission, the scattering of diffuse starlight, and the 
depletion of refractory elements in the ISM compared 
to their expected solar abundances (Savage & Sembach 
1996; Draine 2003; Jenkins 2009). During and after its 
injection into the ISM the newly-formed dust is subject 
to various physical processes, including: grain destruc- 
tion by thermal and kinetic sputtering and by vaporizing 
grain-grain collisions in SN-generated shocks, and grain 
growth by cold accretion. The first models for the evo- 
lution of interstellar dust already identified the discrep- 
ancy between the formation rate of dust in the differ- 
ent stellar sources, and their more rapid destruction rate 
in the ISM (Dwek & Scalo 1979, 1980). Grain growth 
by accretion, or a reduction in the grain destruction 
rates were suggested as possible ways to account for the 
observed abundance of interstellar dust (Dwek & Scalo 
1979, 1980; Dwek 1998; Tielens 1998; Zhukovska et al. 


2008; Calura et al. 2010). 

Further attempts to resolve this discrepancy have led 
to a recent revaluation of the destruction rates of dust 
in the ISM (Jones & Nuth 2011; Bocchio et al. 2014; 
Slavin et al. 2015). The study of Slavin et al. (2015) de- 
rived a silicate lifetime of ~ 2 — 3 Gyr, a ~ 2 —3 fold 
increase compared to previous studies. The discrepancy 
between the formation and destruction rates of silicates 
is therefore still subjected to change. Carbon dust, may 
be more efficiently destroyed, requiring its growth and 
re-formation in the ISM (Jones et al. 2013). 

In this paper we show that iron provides the most ro- 
bust case for the need of the growth of grains by accretion 
in the ISM. We show that because of its unique origin, 
iron is mostly injected into the ISM in gaseous form. 
Yet depletion studies show that most of the iron is ab- 
sent from the ISM gas, and must therefore be locked up 
in dust. This discrepancy is not affected by uncertain- 
ties in the dust destruction rates, and therefore provides 
conclusive proof that most of the iron must have been in- 
corporated into dust by cold accretion in the dense ISM. 

The paper is organized as follows. In Section 2 dis- 
cuss the origin of iron which, unlike other refractory el- 


ements, is predominantly formed in SN Ia. The basic 
equations for following the evolution of the refractory el- 
ements of Mg, Si, and Fe in the solar neighborhood are 
discussed in Section 3. In Section 4 we calculate the max- 
imum amount of Mg, Si, and Fe that can be locked up 
in thermally-condensed grains. Section 5 compares the 
observed depletion trends of these elements, illustrating 
the distinct growth history of Fe bearing grains. The 
results of the paper are briefly summarized in Section 6. 


2. THE ORIGIN OF IRON 


Iron is the most tightly bound nucleus in nature with 
a binding energy of 8.8 MeV per nucleon, and therefore 
the heaviest element that can be synthesized by the fu- 
sion of lighter nuclei in massive stars. It has four sta- 
ble isotopes (A = 54, 56, 57, and 58), with almost all 
(~ 92%) of its observed natural abundance in °°Fe. The 
iron that is injected into the interstellar medium (ISM) 
is made from material that underwent hydrostatic O and 
Si burning. Its isotopic signature is determined by the 
abundance of available free particles (n, p, and a) gener- 
ated in the photodissociation of ?8Si (silicon burning) 
and finalized in the subsequent explosion. It can be 
synthesized either directly as °°Fe or as its radioactive 
precursor °°Ni depending on the neutron excess ratio, 
n = (N — Z)/A, in the expanding material. The produc- 
tion of °°Fe will be dominant when n ~ 4/56 ~ 0.071, 
whereas °°Ni will be the dominant specie at low values of 
n (Clayton 1968; Arnett 1996). The photodissociation of 
8Si occurs at temperatures that are high enough so that 
n = 0, and °6Ni is synthesized by the assembly of a par- 
ticles. The energy released in the radioactive chain °°Ni 
(71/2 = 6.1d) —+ *©Co (11/2 = 77.2d) —> °°Fe provides 
most of the energy that powers Type Ia and CCSN light 
curves. 

Iron is primarily produced in SN Ia events, caused by 
the explosion of a carbon-oxygen white dwarf that ei- 
ther accreted matter from a companion red giant star, 
or that merged with a white dwarf companion. It is 
also produced at the endpoint of the quiescent evo- 
lution of massive stars that end their life as CCSN. 
A 25 Mo supernova (SN) produces only about 0.07- 
0.23 Mo of iron, depending on various parameters such 
as the initial stellar metallicity, the internal stellar ro- 
tation and mixing, and the characteristics of the pis- 
ton (its location in mass and its energy) that simulates 
the explosion (Nomoto et al. 2006; Heger et al. 2003; 
Heger & Woosley 2010; Nomoto et al. 2013). SN Ia pro- 
duce significantly more iron, theoretically about 0.32- 
1.1 Mo, depending on the ignition sites, and the geome- 
try and propagation of the nuclear flame (‘Travaglio et al. 
2004, 2005; Seitenzahl et al. 2013) and observationally 
clustered around 0.4 and 1-1.4 Mg, depending on the 
explosion mechanisms (Childress et al. 2015). 


3. THE CHEMICAL EVOLUTION OF THE SOLAR 
NEIGHBORHOOD 


We use the analytical model described in 
Dwek & Cherchneff (2011) to follow the chemical 
evolution of the solar neighborhood, defined as the 
region contained within a ~ 1 kpc radius around the 
sun. The solar neighborhood has been the subject of 
intensive studies in the past, starting with the pioneer- 
ing work of ‘Tinsley (1974) and more recent studies by 
(Dwek & Scalo 1979, 1980; Dwek 1998; Zhukovska et al. 
2008; Calura et al. 2010, and references therein). Here 
we concentrate on the evolution of Fe, as well as the 
main other, non-carbonaceous, refractory elements Mg 
and Si using the observational constraints on the stellar 
and gas masses, and the star formation, CCSN, and 
SN Ia rates. We adopt the Kroupa stellar initial mass 
function (IMF) in all the calculations (Kroupa 2001). 


3.1. Stellar and Interstellar Gas Mass and the Star 
Formation Rate 


The evolution of the gas mass surface density, U,, in 
the instantaneous recycling approximation is given by: 
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where R-; = 0.50 is the fraction of the initial stellar 


mass that is returned back to the ISM, calculated by av- 
eraging the mass-dependent ejecta over the stellar IMF, 
and where the star formation rate (SFR) surface den- 
sity, Us y(t), is given by the Schmidt-Kennicutt relation 
(Kennicutt & Evans 2012): 


E4(t) fF 
ds r(t) = ds r(t 7 ’ 2 
jolt) = Exsrtto) [SHO (2) 
where Us fr(to) Ug(to) are, respectively, the SFR and gas 
mass surface densities at the current epoch, to = 13 Gyr, 
and where k = 1.5. The metal-free infall rate was taken 
to be an exponential function: 


(Fe), = (=24) exp(-t/r) . (3) 


The stellar mass surface density at time t¢ is simply given 
by the difference: 


x(t) = Ying [1 —ewp(—t/r)|-Xy(t) - (4) 


In Equation (3), Sin ¢ and 7 are model parameters chosen 
to fit the current observational constraints on the surface 
densities of the current SFR, gas, and stellar mass in the 
solar neighborhood. The values of ©;,¢ and 7 used in the 
fits were 40 Mo pe? and 4 Gyr, respectively. Figure 1 
presents the evolution of Ug, Ustar, and Use, in the solar 
neighborhood as as function of time. The observational 


constraints at the current epoch are indicated by vertical 
lines . 
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Figure 1. The evolution of the surface densities of the star for- 
mation rate, and the mass of stars and the ISM in the solar neigh- 
borhood. The vertical bars represent the observational constraints 


on their current values. Details in text. 


3.2. The Evolution of Mg, Si, and Fe 


The Mg, Si, and Fe enrichment of the ISM is domi- 
nated by their synthesis in supernovae, with stellar winds 
primarily returning the mass of preexisting heavy ele- 
ments back to the ISM. The Mg, Si, and Fe masses 
returned from AGB stars were calculated by multiply- 
ing the returned stellar masses, calculated for exam- 
ple by van den Hoek & Groenewegen (1997), by the so- 
lar abundances given by Asplund et al. (2009). CCSNe 
yields, calculated for solar metallicities, were taken from 
Heger & Woosley (2010), and scaled by a factor of 0.60 
to take their dependence on metallicity into account. 

SN Ia yields depend on the location of the detonation 
point in the white dwarf interior. Seitenzahl et al. (2013) 
calculated the nucleosynthesis yield of SNe Ia using 3- 
dimensional models for a range of ignition geometries 
characterized by the density and symmetry of the igni- 
tion points. We adopt here the range of Mg, Si, and Fe 
yields from models N100, N100L and NH100H, charac- 
terized by 100 centralized ignition centers. The nominal 
Fe yield of this model is 0.64 Mo, with low and high val- 
ues of 0.55 and 0.77 Mo. For comparison, an off-center 
model with fewer ignition centers, model N3, produced 
1.1 Mo of iron. Observations of [Co TIT]A5893 emission 
during the nebular phase of the evolution of 7 SNe Ia 
show that the masses of °°Ni synthesized in the ejecta 
cluster around 0.4 and 0.6-1.2 Mo, depending on the 
width of the SN light curve (Childress et al. 2015). The 
range of Fe yields adopted here is therefore consistent 
with the observed range of Fe yields for a mixture of SN 
light curves. 
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An additional factor determining the relative amount 
of Fe synthesized in core collapse and Type Ia SNe is the 
relative rate of these events. The progenitors of SN Ia 
are in the ~ 2— 8 Mo mass range (Maoz 2008; Greggio 
2010), and those of CCSN are stars more massive than 
8 Mo (Smartt 2009; Langer 2012). Consequently, the rel- 
ative contribution of the two sources depends on the stel- 
lar initial mass function (IMF). Because low-mass stars 
have a longer main sequence lifetime than their more 
massive counterparts, the relative contribution of SN Ia 
and CCSN to the Fe abundance depends on the star for- 
mation history as well. Here we adopted the observa- 
tional constraints of the relative SN Ia and CCSN rates 
summarized by Zhukovska et al. (2008) and Calura et al. 
(2010). The nominal ratio of these rates varies from 0.15 
to 0.25. 

Figure 2 (top panel) depicts the evolution of the mass 
of Mg, Si, and Fe, normalized to their solar values, as a 
function of time. A low value of 0.15 was adopted for 
the SN Ia/CCSN ratio of the rates, which reproduces 
the solar abundances at the epoch of the Sun’s birth at 
t = 8.47 Gyr fairly well. The bottom panel shows the 
contribution of the different stellar sources to the total 
abundances of these elements at the current epoch, nor- 
malized to the Fe contribution from SN Ia. The figure 
illustrates the dominance of SN Ia as the source of inter- 
stellar iron. 


4. THE MAXIMUM FRACTION OF MG, SI, AND 
FE IN THERMALLY-CONDENSED DUST 


4.1. Evidence for Dust in CCSNe Ejecta 


SN condensed dust can be observed either during its 
formation process, or later, by looking at the debris of 
the explosion before it has been completely mixed with 
the ISM. The most prominent example for dust for- 
mation in SNe is SN1987A, where ~ 0.4 Mo of dust 
has been inferred from the infrared (IR) emission de- 
tected about 20 yr after the explosion Matsuura et al. 
(2011, 2015). About 0.1-0.2 Mo of SN condensed dust 
has been observed in the remnants of the Crab nebula 
(Gomez et al. 2012; Temim et al. 2012; Temim & Dwek 
2013), and about ~ 0.1 Moe in that of Cassiopeia A 
(Rho et al. 2008; Barlow et al. 2010; Arendt et al. 2014). 
A detailed review of searches for the evidence of dust for- 
mation in SNe has been presented by Gall ct al. (2011). 

The observed dust mass in CCSNe falls always short 
of the total mass of refractory elements present in their 
ejecta. In particular, most of the iron produced in CCSN 
is observed to be in the gas phase. A survey of the 
X-ray emitting ejecta of Cas A gives an Fe mass of 
~ 0.09—0.13 Mo in the hot gas (Hwang & Laming 2012), 
comparable to the total mass of iron synthesized in the 
explosion. This value depends on the assumed volume 


filling factor of the emitting region, and should therefore 
be regarded as an upper limit. In any case, the very pres- 
ence of iron in the X-ray emitting gas shows that not all 
of it condensed in the ejecta. In spite of these observa- 
tions, we will assume that all the Fe produced in CCSNe 
condenses into dust. 


4.2. Evidence for the Absence of Dust in SN Ia Ejecta 


Type Ia SN are less likely to form dust in their ejecta 
compared to CCSNe. The mass ejected in a SN Ia 
event is typically ~ 1.4 Mo, significantly lower that the 
~ 5-15 Mo of material ejected in a typical CCSN. With 
similar explosion energies and expansion velocities, the 
temperature and density of the SN Ia ejecta decreases 
more rapidly compared to those in CCSNe ejecta, short- 
ening the time available for the formation of molecules 
and grains. Furthermore, the abundance and growth of 
the grains in SN Ia ejecta is greatly impeded by the large 
abundance of radioactive °°Ni, which produces fast elec- 
trons and y-rays that are effective at dissociating precur- 
sor molecules in this environment. Models for the forma- 
tion of dust in SN Ia show that any dust that forms will 
have small radii (< 100 A), and will probably be totally 
destroyed by the reverse shock which is generated as the 
ejecta is decelerated by the ambient ISM (Nozawa et al. 
2011). The models produced no significant amount of 
iron dust. 

Extensive searches for newly-formed dust in the rem- 
nants of SN Ia have been conducted with the Spitzer 
and Herschel satellites. Observations of remnants 
such as Kepler, (Blair et al. 2007; Williams et al. 2012; 
Gomez et al. 2012), RCW 86 (Williams et al. 2011), 
SN 1006 (Winkler et al. 2013), Tycho (Williams et al. 
2013; Gomez et al. 2012), and N103B (Williams et al. 
2014) have shown that all the IR emission from these 
remnants arises from shocked interstellar dust, with no 
evidence for any dust associated with the SN ejecta 
(Gomez et al. 2012). A firm upper limit that is consid- 
erably less than 0.1 Mo was put on any dust present in 
the ejecta of Kepler’s SNR (Blair et al. 2007). For all 
practical purposes, we can safely conclude that all the 
Fe produced in SN Ia is returned to the ISM in gaseous 
form. 


4.3. Dust in Stellar Winds 


AGB stars are important sources of dust, forming car- 
bonaceous dust when the C/O ratio in their ejecta ex- 
ceeds unity, and forming silicate dust otherwise (Waters 
2004). For AGB stars with initial solar metallicity, sil- 
icates will form in stars with masses in the 1.0-1.7 and 
5.0-8.0 Meo mass range (Nanni et al. 2013). Here we as- 
sume that all the Mg, Si, and Fe in the ejecta of stars in 
this mass range are locked up in dust. 
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Figure 2. Top panel: The evolution of Mg, Si, and Fe abun- 
dances, normalized to solar abundances, versus time. The horizon- 
tal dashed lines indicate deviations by a factor of two from solar 
abundances (thick horizontal line). The hatched area represents 
the range of nucleosynthetic yields of these elements in SN Ia. The 
vertical solid line at 8.47 Gyr marks the age of the sun. Bottom 
panel: The relative contribution of the different sources (AGB 
stars, CCSNe, and SN Ia) to the total abundance of Mg, Si, and 
Fe, at the current epoch (13 Gyr), normalized to the SN Ia contri- 
bution to the Fe abundance. The error bars represent the adopted 


range of Fe yields in SN Ia. 


4.4. Maximal Fraction of Mg, Si, and Fe in 
Thermally-Condensed Dust 


In calculating the maximum possible depletion of Mg, 
Si, and Fe in thermally condensed dust, we will assume 
that they condense with unit efficiency in the ejecta of 
CCSN and AGB stars, but that they are returned in 
gaseous form in the ejecta of SN Ia. This will provide 
a strong upper limit on the fraction of Fe that can be 
locked up in thermally condensed dust in the ISM. 

The results of our calculations are presented in Fig- 
ure 3 which depicts the value of faust(X ), defined as the 
fraction of the abundance of the element X ={Mg, Si, 
Fe} that is locked up in thermally-condensed dust, 


faust (X) = Xaust/Xtot ’ (5) 


where X gust is the abundance of X in the dust, and Xyoz 
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Figure 3. The calculated fraction, fqgust(X ) of the total abun- 
dance of X = {Mg, Si, Fe} that is locked up in thermally-condensed 
dust (grey bars), is compared to observed values in the cool (blue) 
and hot (red) phases of the interstellar medium (Savage & Sembach 
1996). 


is its total (gas+dust) abundance. A value of faust(X) = 
1 indicates that the element X is totally absent from the 
gas phase of the ISM. The calculated value of faust(X) 
represents the maximum attainable depletions of an el- 
ement by thermal condensation, since it assumes a con- 
densation efficiency of 100% in AGB stars and CCSNe, 
and that none of the dust formed in CCSN is destroyed 
by the reverse shock generated by the propagation of 
the SN blast wave into its surrounding medium (Dwek 
2005; Bianchi & Schneider 2007; Nozawa et al. 2007; 
Nath et al. 2008; Biscaro & Cherchneff 2016; Silvia et al. 
2010, 2012; Micelotta et al. 2016). 

Comparison of the calculated values of faus:(X) to the 
observed depletion fractions in the cool and warm phases 
of the ISM (Savage & Sembach 1996) shows that, in gen- 
eral, they fall short of their depletion fractions in either of 
these phases of the ISM. The differences are not so large 
for Mg and Si bearing grains. The lifetime of silicate dust 
is still a subject of intense research, and may further in- 
crease when the clumpy nature of the ISM is taken into 
account (Dwek et al. 2007; Jones & Nuth 2011). The ob- 
served depletions of Mg and Si may yet be explained by 
the efficiency of dust formation in the sources. The dif- 
ference between the maximum attainable depletion and 
the observed depletions in the ISM is most pronounced 
for Fe. At most, only 35% of the iron can be locked up in 
dust, compared to the observed 99 and 94% depletions 
in, respectively, the cool and warm phases of the ISM. 

X-ray absorption line spectra provide an additional 
means for distinguishing between the gas and solid phase 
abundances of select refractory elements, and the com- 
position of dust in the ISM (Lee et al. 2009). Indeed, 
X-ray absorption line studies towards eight X-ray bina- 
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Figure 4. Results of absorption line studies in the soft X-ray 


band region with the Chandra spacecraft (dots). The green solid 
line is a parametric fit using metallic iron data. The red line repre- 
sents the iron absorption edges if all the iron was in the gas phase. 
Data and the parametric fit were provided by T. Kallman, private 


communication. 


ries, obtained with the Chandra High Energy Transmis- 
sion Grating Specrometer provided evidence for the fact 
that most of the interstellar iron is locked up in dust 
(Gatuzz et al. 2015). Figure 4 shows the observed soft X- 
ray spectrum around the location of the iron L2 and L3 
absorption edges. The observations are represented by 
dots, and the green line represents a model in which the 
cross section and absorption edges were obtained from 
metallic iron data. The data could not be fitted with a 
pure iron gas model which is characterized by sharp ab- 
sorption edges as shown by the scaled red line. The X-ray 
data provide therefore additional evidence that most of 
the interstellar iron is locked up in dust. So even if the Fe- 
bearing CCSN and AGB dust grains are not destroyed in 
the ISM, the very low faust(F'e) value compared to the 
observed depletions provides the most compelling evi- 
dence for the need for its depletion in the ISM. 


5. MAGNESIUM, SILICON, AND IRON 
DEPLETION TRENDS 


Additional insight into the nature of Fe-bearing dust 
can be obtained from a comparison between its deple- 
tion to that of other, non-carbonaceous, refractory ele- 
ments. Carbon dust follows a distinct evolutionary cycle 
(Jones et al. 2013), and the analysis of its depletion pat- 
tern is outside the scope of this paper. The depletion of 
an element X from the gas phase is usually defined as: 
[X (gas)/H] = log{N(X)/N(H)} — log(X/H)o, which 
expresses the reduction in its gas phase abundance (rel- 
ative to hydrogen) compared to the expected value if all 
the X atoms were in the gas phase. N(X) and N(H) 
are, respectively the observed abundance of the element 
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Figure 5. Top panel Interstellar depletion factors [Fe(gas) /H] 
versus [Mg(gas)/H] and [Si(gas)/H]. Bottom panel The fraction 
of the total abundance of Mg and Si locked up in dust versus the 
Fe fraction locked up in dust. Curves were calculated using the 
elemental depletion parameters of Jenkins (2009). The dashed line 
shows the ISM depletion trend if Fe was depleted at the same rate 


or with the same stoichiometric ratio as Mg and Si. 


X and H (in atomic plus molecular form) along the LOS. 

A detailed study of the depletion pattern of 17 different 
elements was conducted by Jenkins (2009). Surveying 
over 243 lines of sight (LOS), he defined a parameter F* 
to represents the overall strength of the depletion along 
any given LOS. A value of F* = 0 was assigned to a 
LOS with the lowest overall depletion, whereas a value of 
F* = 1 was assigned to the LOS toward the star ¢ Oph, 
known for its strong depletions. Though most F™* values 
fall in the {0,1} interval, some LOS are characterized 
by lower or larger values. The value of F™* reflects the 
history of grain growth and destruction in the gas along 
a given LOS. 

Figure 5 (top panel) plots the observed depletion of 
Fe versus that of Mg and Si. The bottom panel de- 
picts the same information in the form of the fraction 


of the different elements locked up in the dust, given 
by (X(dust)/H)/(X/H)o = [1 — 10!*(925)/4]], The el- 
ements Mg, Si, and Fe attain the same maximum ob- 
served depletion, [X(gas)/H], for F* = 1 in the ISM, 
with about 95% of the Mg and Si, and 99% of the Fe 
locked up in dust. However, the minimal observed de- 
pletion, [X(gas)/H]o for F* = 0, of Mg and Si is about 
~ 40 — 45%, whereas that of Fe is ~ 89%. This means 
that any LOS in the ISM represents material in which 
the depletion of Fe has already grown from ~ 35 in the 
sources to ~ 90% in the ISM. Also shown in Figure 5 
(bottom panel) is the correlation between the depletions 
of Mg, Si, and Fe if they were depleted to form minerals 
of the olivine group with a stochiometric ratio of 1:1:1 
(dashed curve). The observed trends (colored lines) sug- 
gest that the progress of the depletion of Mg and Si is 
unrelated to that of Fe. If Mg, Si, and Fe accreted onto 
grains at this stochiometric ratio, then only ~ 50% of the 
Mg and Si abundance could be locked up in dust when all 
the gas phase iron is exhausted. Mg and Si must there- 
fore follow a distinct accretion process. Furthermore, the 
significantly larger observed value of |Fe(gas)/H]o com- 
pared to that of Mg and Si, in spite of its initial low value 
of fin(Fe), shows that Fe accretes more efficiently onto 
the pre-existing thermally-condensed grains than either 
Mg or Si. 

Support for the different depletion trends is provided 
by recent spectral analysis of the composition of inter- 
stellar grains towards ¢ Ophiuchi, a dense interstellar 
cloud representing the transition between the diffuse and 
dense phases of the ISM (Poteet et al. 2015). The ob- 
servations show that Mg and Si are depleted in silicate 
compounds, whereas Fe is depleted in distinctly differ- 
ent solids. These observations provide important con- 
straint on any model for an ISM origin for interstellar 
dust, namely preventing the accretion of iron onto sili- 
cate compounds. 


6. SUMMARY 


Recent re-evaluation of dust lifetimes in the ISM con- 
cluded that the destruction rate of silicate dust may 
have been overestimated in the past, significantly nar- 
rowing the discrepancy between their formation rates in 
the different stellar sources, and their destruction rates 
in the ISM (Jones & Nuth 2011; Bocchio et al. 2014; 
Slavin et al. 2015). The need to reconstitute silicates by 
accretion in the ISM is therefore still a subject of further 
studies. 

In this paper we have shown that most of the iron 
depletion must take place in the ISM. This conclusion 
is robust, arising from the fact that most of the Fe is 
formed in SN Ia, and therefore ejected into the ISM in 
gaseous form. The maximum percentage of Fe that can 
be depleted onto thermally-condensed dust is less 35%, 


since we assumed a maximal condensation efficiency in 
AGB stars and CCSNe. 

The general issue of grain growth in the ISM invites a 
new set of well known problems that have to be reconciled 
with observational constraints (Jones & Nuth 2011). An 
additional problem is the ubiquitous presence of the 
9.7 and 18 wm silicate absorption features towards the 
Galactic Center (McCarthy et al. 1980; Fritz et al. 2011; 
Roche et al. 2015, and references therein). The features 
are attributed to, respectively, the Si-O stretching and O- 
Si-O bending modes in silicates comprised of SiO, tetra- 
hedra (Khanna et al. 1981). Forming such tetrahedral 
structures requires the silicates to have formed or an- 
nealed at temperatures above ~ 1000 K for tens of hours 
(Hallenbeck et al. 2000). These conditions are unattain- 
able if most of the silicates are grown by cold accretion 
in the ISM. Interstellar dust must therefore mostly com- 
prise of composite refractory grains, carbon and silicates 
with organic refractory inclusions (Li & Draine 2001; 
Zubko et al. 2004), or composite silicate carbon grains 
(Jones et al. 2013). Reconciling these requirements with 
the observed discrepancies between dust formation and 
destruction rates in the local ISM and galaxies remains 
a major problem in astrophysics. 

Our work and observations by Jenkins (2009) and 
Poteet et al. (2015) suggest that grains with accreted 
iron inclusions should be considered as a separate dust 
component of the ISM. Recent in-situ analysis of in- 
terstellar dust that penetrated the solar system sup- 
ports this picture. The Cosmic Dust Analyzer on board 
the Cassini spacecraft revealed that the passing inter- 
stellar dust grains consist primarily of magnesium-rich 
grains of silicate oxide composition with iron inclusions 
(Altobelli et al. 2016). 

Tron inclusion in grains will have important effects 
on their alignment properties (Mathis 1986). Labora- 
tory measurements of the optical properties of coated 
and composite dust particles, and exposure of compos- 
ite grains containing metallic and organic material to 
UV and particle irradiation are important for furthering 
our understanding of the origin and nature of interstellar 
dust. 


In writing this paper I have benefitted from many en- 
lightening conversation with Joe Nuth, and useful refer- 
ences provided by Steven Rodney. J.N. and Rick Arendt 
provided useful comments on an early version of the 
manuscript. I thank the referees Xander Tielens and An- 
thony Jones for their critical comments which led to im- 
provements in the manuscript. This work was supported 
by NASA’s 12-ADP12-0145 and 13-ADAP13-0094 re- 
search grants. 
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